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Do high-velocity clouds form by thermal instability? 
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ABSTRACT 

We examine the proposal that the H i "high- velocity" clouds (hvcs) surrounding the 
Milky Way and other disc galaxies form by condensation of the hot galactic corona via 
thermal instability. Under the assumption that the galactic corona is well represented 
by a non-rotating, stratified atmosphere, we find that for this formation mechanism to 
work the corona must have an almost perfectly fiat entropy profile. In all other cases 
the growth of thermal perturbations is suppressed by a combination of buoyancy and 
thermal conduction. Even if the entropy profile were nearly flat, cold clouds with sizes 
smaller than 10 kpc could form in the corona of the Milky Way only at radii larger than 
100 kpc, in contradiction with the determined distances of the largest HVC complexes. 
Clouds with sizes of a few kpc can form in the inner halo only in low-mass systems. 
We conclude that unless even slow rotation qualitatively changes the dynamics of a 
corona, thermal instability is unlikely to be a viable mechanism for formation of cold 
clouds around disc galaxies. 

Key words: ISM: kinematics and dynamics, galaxies: haloes, galaxies: kinematics 
and dynamics, galaxies: evolution 



1 INTRODUCTION 

It is widely accepted that star-forming disc galaxies such 
as the Milky Way have sustained their star formation 
over gigayears by accre ting gas at a fairly steady or 
slowl y declining rate (e.g. IChiappini. Matteucci fc Romanq 
I2OOII ). The source of this gas is uncertain. The Milky 
Way and other star-forming disc galaxies (such as M31) 
that have been studied with sufficient sensitivity in 



servation of CIV, OV and OVI absorption along lines of 
sight that pass close to the Galax y's HVCs (|Sembach et al.l 
l2003l : lFox. Savage fc Wakkeill2006l ). and the morphology of 
the Magellanic Stream and individual HVCs jPutman et al.l 
l2003l : iBriins et al.ll2000l ) . 

The relationship of the whim to circum-galactic clouds 
of Hi is unclear. Undoubtedly much of the Hi that is ob- 
served more than a kiloparsec from galactic planes has 
been pushed off the star-forming pla ne by the mechani- 



rounded by clouds of H i (e.g. 


Wakker & van Woerdenll 19971: 


IWestmeier, Braun & Thilkerl 


20051). but the amount of gas 



contained in these "high-velocity" clouds (hvcs) is sufficient 
to sustain accretion onto the disc for ab out a gigayear (e.g. 
IWakker et all l2007l : ISancisi et a"l]|2008l ). The only realistic 
source of gas for sustained accretion onto discs is the warm- 
hot intergalactic medium (whim), which is believed to con- 
tain over half the baryons in the Universe. This belief is 
based on a combination of big-bang nucleosynthesis theory 
and observations of the cosmic m icrowave background (e.g. 
iFukugita. Hogan fc Peebleslll998l ). so within the framework 
of standard cosmology it must be considered robust. Em- 
pirical evidence for the existence of the whim includes the 
requirement for a medium to con fine clouds of i nterstellar 
gas far from the Galactic plane (|Spitzer I Il956l ) , the ob- 
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cal fe edback from star formation (e.g. Fraternali fc BinnevI 



2008h - t h e gala ctic fountain effect (|Shapiro fc Fieldlll976l ') 



BregmanI (|l98(]| ') proposed that the galactic fountain model 



could explain also the formation of the H VCs by condensation 
from hot gas outflows from the disk, but iFerrara fc Einaudil 
( 1992J argued on theoretical grounds that this mechanism is 
unlikely to work. The hypothesis that HVCs formed by con- 
densation of the hot gas of the galactic fountain is discarded 
also on observational grounds, because such gas cannot be 
less metal rich than the local interstellar medium, whereas 
the high- velocity cloud Complex C, which lies 6 — 11 kpc 
from the plane, has Z ~ 0.15Zq l|Wakker et al.ll2007^ . This 
finding suggests that some HVCs, especially those more than 
~ 5 kpc from the plane, repr esent matt er that is entering 
the Galaxy for the first time l|Oortlll970l) . 

At least some HVCs must have been stripped from in- 
falling satellites because some clouds are clearly associ- 
ated with the Magellanic Stream, which is itself thought 
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to have been torn from the Small Magellanic Cloud 
l|Putman et aL 20031). Count s of stars in the Sloan Digital 
Sky Survey ^elT et alj lioosl ) and in deep surveys of M31 
l|Richardson et al.. .2008 ) have shown that the outer stellar 
halos of the Galaxy and M31 are far from relaxed and are 
best thought of as a superpositions of dissolving satellites. 
Any gas that belonged to these satellites is likely to give rise 
to HVCs. 

However, the supply of cold gas from infalling satellites 
falls far short of what a typical star- forming disc galax y 
needs to sustain its star formation (e.g. ISancisi et al.ll2008l ). 
In the long run that gas must come from the whim because 
the latter is the only gas reservoir of sufficient capacity. So if 
star-forming galaxies obtain the gas they need by accreting 
HVCs, these clouds must form from the whim. A case for the 
formation of HVCs from the whim has been argued on the ba- 
sis of simula tions of galaxy formation that follow the baryons 
with SPH llKaufmann et al.1 [ioO^: ISommer-LarsenI [200^: 



Peek. Putman fc Sommer-Larsen 120081 : Kaufmann et alJ 



20091 ). This case is not strong because (a) the mass scale of 



HVCs is smaller than the mass resolution of state of the art 
cosmological simulations, (b) there are grave doubts about 
the reliability of SPH when the fluid contains extremely 
steep density gradients such as those expected at the inter- 
face of the WHIM and an H I cloud and (c) multi-phase media 
in SPH sim ulation s might be subject to artificial overcooling 
dMarri fc White 200^ : iKaufmann et al.1 l2007l : lAgertz et at] 
I2OO7I . but see also Price 2008). In this paper we use largely 
analytic arguments to examine the likelihood that H I clouds 
can condense from the whim > 10 kpc above the main H I 
disc. 

The dynamics of the whim around a galaxy such as the 
Milky Way bears a close resemblance to the dynamics of 
"cooling flows" in groups and clusters of galaxies. The gas 
in a cooling flow is hotter and denser than the whim, so 
it can be readily studied through its X-ray emission. When 
observations showed that the cooling time at the centre of 
a typical cooling flow is significantly less than the Hubble 
time, it was assumed that the gas was flowing towards the 
centre and th ere generating a pool of cold gas from which 
stars formed l|Silljil976l : [Cowie fc Binnevll 19771 ). Shortly af- 
terwards it became clear that this picture was not viable 
because (a) it predicted a radial X-ray surface-brightness 
profile that was too centrally concentrated, and (b) insuffi- 
cient young stars were present at small radii, no matter what 
the initial mass function. The hypothesis then took firm root 
that ~ 100 kpc from the centre, thermal instability i n the in- 
flowi ng gas caused clouds of cold gas to "drop out" (|Nulsenl 
Il986l ). These clouds were presumed to be too small to be de- 
tected, but removed sufficient X-ray emitting gas from the 
cooling flow to reconcile the theoretical and observational 
surface-brightness proflles. Eventually this theory of "dis- 
tributed mass drop out" was ruled out by the absence of 
spectral lines in soft X-rays that should have ac companied 
the formation of the clouds (|Peterson et al.ir2002l ). However, 
more than a decade before the decisive observational evi- 
dence arrived the in tellectual foundations of the theo ry had 
been shot away by iMalagoli. Rosner fc Bodol l| 19871 . here- 
after MRB), who pointed out that a hot atmosphere that 
is confined by a gravitational field and radially stratified 
by specific ent ropy is not thermally unstable in the sense 
of iFieldl (|l965l ). Subsequent papers clarified important as- 



pects of the problem (Balb us fc Soker 1989, hereafter BS, 
lTribble|[l989l : iBalbud [19911 ) but confirmed the basic phys- 
ical principle. Unfortunately, it seems that the community 
that is now w orking on galaxy forma tion is unaware of this 
principle fe.g. lMaller fc BuUocklliooi . Peek et al. 2008), and 
there is a danger that studies of the whim will fall into the 
same trap as did early studies of cooling flows. In particular, 
condensation of the whim into HVCs > 10 kpc above the H I 
disc is an instance of distributed mass drop out, albeit in a 
different parameter regime. In this paper we examine this 
possibility critically. 

In Section 2 we summarise the analytical work of MRB 
and BS, and explain the underlying physics. In Section 3 we 
describe the simple models of galactic coronae in which we 
study thermal instability. Section 4 presents both the result 
of applying linear theory to these models and studying the 
time evolution of spherical model atmospheres. Section 5 
sums up and discusses the implications of our results for 
galaxy accretion. 



2 THERMAL INSTABILITY IN STRATIFIED 
CORONAE 

The WHIM around a galaxy like the Milky Way is expected 
to be approximately in equilibrium, stratified in the galactic 
gravitational potential. The thermal instability in gravita- 
tionally stratifie d coronae ha s been stu died by several au- 
thors (e.g. MRB. IWhite fc Sarazin.l987l : lTribbl6jil989„ . BS), 
mainly in the context of the study of cooling fiows in galaxy 
clusters. In many respects a galactic corona is just a scaled- 
down version of the hot atmosphere of a galaxy cluster. 
Thus, under the assumption that the corona is non-rotating 
and spherically symmetric, the formalism of previous work 
of spherical models of cooling flows can be applied straight- 
forwardly. However, the detailed thermal behaviour of the 
system is sensitive to the gas temperature and density dis- 
tribution, so the conclusions drawn in the case of clusters 
cannot be a priori extended to lower-mass systems. 

Here we follow the treatment of MRB who used Eule- 
rian plane-wave perturbations to study thermal instability 
in a spherically symmetric stratified atmosphere in the pres- 
ence of a fiow, including the effects of cooling and thermal 
conduction. It has been pointed out that when a background 
ffow is present, the use of Lagrangian p erturbations is prefer- 
able to that of Eulerian perturbations (Balbus 1988;'Tribbl^ 
il989i , BS) . However, the main results of MRB have been con- 
firmed by the Lagrangian study of BS. Thus, for the purpose 
of the present investigation we prefer to use the Eulerian ap- 
proach, which has the advantage of a simpler formalism than 
the Lagrangian approach. 



2.1 Perturbation analysis 

We summarise here the perturbation analysis of a non- 
rotating spherical stratified atmosphere by MRB. The sys- 
tem is governed by the equations for mass, momentum and 
energy conservation 

dp 
dt 

9v 



+ V • (pv) = 0, 



— +v Vv = ^ 

at p 



(1) 
(2) 
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A(r). 



(3) 



Here $ is the galactic gravitational potential, 7 = 5/3 is 
the ratio of principal specific heats, A is the cooling func- 
tion, jj, is the mean gas particle mass in units of the proton 
mass TTip, Ks — 1.84 x 10~^ ergs~^ cm~^ K"'^''^ is Spitzer's 
coefficient of thermal co nductivity (as suming a value 30 for 
the Coulomb logarithm; ISpitzeij|l962l l and / < 1 is the fac- 
tor by which therm al conduction is suppre ssed by a tangled 
magnetic field (e.g. iBinnev fc Cowid [l98ll ) . We are assum- 
ing that the thermal conductivity is unsaturated, which will 
certainly be the case in an isothermal corona since there 
the only temperature gradients are those associated with 
perturbations and therefore of order 5T/X, where ST is in- 
finitesimal and A is the perturbation's wavelength. Near the 
edge of an adiabatic corona the temperature gradient is steep 
and the density low (Fig. [S] below), so the conductivity may 
saturate. Reduced conductivity will underline the tendency 
of these regions to be thermally unstable. 

Apart from the presence of the factor /, we treat the 
hot galactic gas as unmagnetized. We do not attempt to 
model the magnetic field in the corona for simplicity and 
because its properties are poorly constrained. An ordered 
ma gnetic field can have subtle effec ts on thermal instabil- 
ity l|Loewenstein|[l990l : lBalbu^ll99lf ). but, as we explain in 
Section 5, far from the disc the field is unlikely to be ordered. 

The unperturbed corona is assumed to be close to hy- 
drostatic and thermal equilibrium in the sense that the sys- 
tem is approximately in a steady state over the time scales 
of interest even in the presence of radiative cooling and ther- 
mal conduction. Thus, the corona can be described by the 
time- independent spherically-symmetric pressure p, density 
p, temperature T and velocity field v — v&r, which satisfy 
equations ((T} to ((31) with vanishing partial derivatives with 
respect to t. 

We linearise equations (|H3|I with Eulerian perturba- 
tions of the form F + (5F exp(— ici^t + ik^r + ikt ■ Xt), where 
xt is the tangential vector, is the radial wave- number, 
kt = |kt| is the tangential wave-number and \5F\ ^ \F\. Un- 
der the assumption of short wave-length and low-frequency 
perturbations, the linearised equations reduce to the follow- 
ing dispersion relation: 



+ 1(D(lJc + tOth) — ^i^BV = 0, 



(4) 



where uj = cj — kiV is the frequency experienced by a pertur- 
bation that rides with the background fiow and k^ = k'^ + k^. 
In the equation above 

7/2 



_ 7 - 1 k^fKsT 

(^c = , 

7 p 

7-1 nA(r)A(r) 

"^th = ; — tt; 1 

7 ksT 

where n — p/ finip is the particle number density. 



A(r) 



din A 
dlnT' 



(5) 
(6) 

(7) 




Figure 1. Contours of constant specific entropy (in the plane 
of the tangential and radial coordinates) for a model corona with 
perturbations that have fct/fcr = 2/7 (upper panel) and 7/2 (lower 
panel). The gravitational field is directed from top to bottom. 
Specific entropy is in arbitrary units. 
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dT /dr\ 
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dT ^ 7-1 g/imp 
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1/2 



1/2 



(8) 



and ljbv is Brunt- Vaisala (hereafter BV) frequency defined 



where g — ||V$|| and (dT/dr)adiab is the adiabatic temper- 
ature gradient. The BV frequency ujbv is clearly real in the 
relevant case of convectively stable configurations. 

In the absence of cooling and thermal conduction (wc = 
i^th ~ 0), the dispersion relation yields 

£l> = ■^^BV. (9) 

Thus ljbv is the frequency at which radially-oriented needles 
of gas bob up and down around their equilibrium radii. As 
kt/k falls from unity, the shape of the oscillating region flat- 
tens through a sphere, reaching a tangentially oriented disc 
in the limit fct/fc — > (Fig.[l|. According to equation ((9|, the 
oscillation frequency declines from ljbv to zero as the body 
is flattened from a needle to a disc. To understand this re- 
sult physically, consider the upper panel of Fig. [1] showing 
isoentropy contours for a model corona with perturbations 
that have kt/ki — 2/7. In the absence of cooling or con- 
duction, speciflc entropy is unchanged as the fluid moves, 
so the colours and contours are frozen into the fluid. At the 
points marked A and C in the flgure, the entropy perturba- 
tion is positive because the fluid is displaced downward from 
equilibrium, while at point B the fluid is displaced upwards. 
The upward movement of the contour above B has clearly 
been enabled by fluid moving away horizontally towards the 
points A and C. In the lower panel - showing isoentropy 
contours for the same model corona, but with kt/kr = 7/2 - 
it is clear that the upward movement at B requires a less sig- 
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nificant horizontal movement of fluid. Consequently, the ef- 
fective inertia that has to be overcome by the perturbation's 
buoyancy increases with the flattening of the perturbation, 
so uj decreases with kt/k. 

In the absence of conduction or a conflning gravitational 
force [loc ~ i^bv ~ 0) equation Q stat es that u) = — icJth, 
so the flow is unstable when t^th < (|Fieldlll965l l. From 
equation ([6} ujth < when A(r) > 0, which is always the 
case at temperatures > 7x 10'' K, as is apparent from Fig.[5J 
in which A(T) is plotted for th e tabulated cooling function 
of [Sutherland fc Dopital l| 19931 ). which we adopt through- 
out the paper. Thus, the condition for thermal instability 
is satisfied in the temperature range of interest and we can 
specialise to the case tJth < 0. The conduction frequency 
ojc > always, and conduction replaces the Field instability 
with exponential decay when ujc > \uJth\- Since loc cc k'^, the 
shortest-wavelength modes are always damped. 

In the course of buoyant oscillations, a parcel of gas is 
half the time overdense and cooler with respect to its sur- 
roundings, and half the time underdense and warmer, so 
unless \u}c + Wthl ^ (kt / k)ijjBv , neither conduction nor radi- 
ation will have a big impact on buoyant oscillations. This 
expectation is borne out by exact analysis of the disper- 
sion relation (|4]), which shows that when ljc = modes are 
monotonically growing if kt/k is smaller than 



2ujbv 



(10) 



In addition, thermal conduction stabilises perturbations 
with wavelengths smaller than 



A ■ - 2^ 

n 



A(r)A(r) 



1/2 



(11) 



Although sufficiently fiattened perturbations can be 
thermally unstable, equation (|lip places a lower bound on 
the thickness of any thermally unstable disc, so such a disc's 
transverse size must be considerable. It is likely that the hot 
atmospheres of galaxies, galaxy groups and clusters are tur- 
bulent by virtue of cosmic infall, galactic winds and jets 
from AGN. Consequently, an extended disc-like overdensity 
is unlikely to remain flat and exactly tangentially oriented 
for more than a dynamical time. Once turbulence has dis- 
torted and/or rotated a perturbation, its value of uj will 
rise towards lobv and it will cease to be thermally unstable. 
Therefore we can discount the possibility that uj <^ ujbv- 



2.2 Physical implications 

The thermal and dynamical evolution of a galactic corona 
is ultimately determined by four timescales: the dynamical 
time of the system t^yn, the gas isobaric cooling time 



^rool — 



7 fceT 
7-1 nA{T) ' 



(12) 



the BV time tev = '-Jbv a-nd the thermal instability 
time tth = I'^^th^l- Figure [3] illustrates this by plotting 
these timescales for a particular atmosphere representa- 
tive of those of massi ve elliptical galaxies (model HM of 
iNipoti fc Binnevll2007l . with metallicity Z = 0.3Zq; see also 
Table [l]). Clearly, a necessary condition for a corona to be 
in equilibrium is tdyn < icooi, so in Fig. |3]the dotted line 




Figure 2. The quantity A(T) = 
of temperature (full li ne: Z = Zp,, 



2 — dlnA/dlnT as a function 
dotted curve Z = O.OSZq) for 



the cooling function of [Sutherland &: Dopital lll993l) 



lies above the short-dashed line: a plasma in which this con- 
dition is violated is experiencing a cooling catastrophe, and 
the question of its thermal stability is meaningless. Let us fo- 
cus on systems with tdyn < icooi , which are approximately in 
equilibrium and, at least in principle, can be prone to ther- 
mal instability. As we have seen, thermal perturbations of 
average shape {kt ~ ki) can grow only if tth < isv, because 
otherwise buoyancy is effective in restoring the unperturbed 
entropy profile. In the typical case of subadiabatic temper- 
ature profile, the BV time is of the order of the dynamical 
time (both being determined mainly by the gravitational 
field), while the thermal instability time is of the order of 
the cooling time (both being determined mainly by the cool- 
ing function). Hence in Fig. [3] the long-dashed line lies close 
to the short-dashed line, while the solid line lies close to the 
dotted line, which necessarily lies above the short-dashed 
line. As a consequence, it is not easy to find systems in 
which both the conditions tdyn < tcooi and tth < tev are 
satisfied, and this is the main reason why stratified atmo- 
spheres are typically thermally stable as MRB pointed out 
in the context of rich galaxy clusters. 

However, from equations (|6]l and (|12p . 

^cool 
tth 



lA(r)|, 



(13) 



and Fig. [2] shows that tth < icooi for r>8 x lO* K. Con- 
sequently, for certain temperature and density distributions 
within the system, tsv can be longer than tth, and tdyn 
shorter than tcooi, even when tev ~ tdyn. Figure [2] shows 
that the scope for simultaneously satisfying both conditions 
is largest at temperatures ~ 3 x 10^ K where the cooling 
function increases steeply for decreasing temperature. Thus, 
the development of thermal instability in galactic coronae, 
which have virial temperatures in this range, cannot be ex- 
cluded a priori. 
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Figure 3. The radial dependence of dynamical, cooling, ther- 
mal instability and BV timescales in a model atmosphere rep- 
resentative of those of massive elliptical galaxies (model HM of 
iNipoti fc Binnevll2007l. see also Table [ij. 



As the temperature profile approaches the adiabatic 
temperature profile, tsv becomes arbitrarily long (see equa- 
tion |8]) - in an adiabatic atmosphere there is no tendency 
for buoyancy to restore displaced material to its original 
radius. Hence an atmosphere with a very shallow entropy 
profile may well be thermally unstable. Information on the 
entropy profile of galactic coronae is then crucial for the 
issue of their thermal stability (see also BS). 

Finally, thermal conductivity can potentially give an 
important contribution to the suppression of small-scale 
perturbations: conductivity damps all perturbations with 
wavelength smaller than the critical wavelength Acrit (equa- 
tion lll[) . The extent of this damping depends on the 
poorly constrained suppression of conductivity, which we 
parametrise with the dimensionless factor /. However, we 
will see that even strongly suppressed thermal conduction 
can have non-negligible stabilising effect. 

The discussion above leads to the conclusion that the 
window for the development of thermal instability in galac- 
tic coronae is narrow, but not necessarily closed. Thus, to 
address the question of the thermal stability of the whim it 
is necessary to apply the results of the perturbation analysis 
to specific models of galactic coronae. 



3 APPLICATIONS TO GALACTIC CORONAE 

Application of the thermal instability analysis to the X-ray 
emitting hot atmospheres of massive elliptical galaxies and 
clusters of galaxies has shown clearly that these systems 
cannot develop thermal instabilities (MRB, BS). This result 
is apparent from Fig. [S] which shows the relevant timescales 
for a massive elliptical galaxy model: in this case the BV 



timescale is at least two orders of magnitude shorter than 
the thermal instability timescale. 

Here we address directly the question of whether the 
WHIM surrounding disc galaxies is prone to thermal insta- 
bility. Unfortunately, the properties of the coronae of disc 
galaxies are not as well known as those of their more massive 
counterparts, mainly because they have low X-ray surface 
brightnesses on account of their low gas densities. Despite 
various attempts to detect in X-rays the coronae of disc 
galaxies, so far only upper limits on their X-ray l uminos- 
ity have been obtained (e.g. iRasmussen et al.|[2009h . In the 
cases in which extra-planar X-ray emission is detected, it ap- 
pears clumpy, with clear connections with the star-forming 
regio ns, and is likely to be explained by g alactic fountains 
(e.g. [Strickland et al.l l2004l : iLi et all [20061 ). But also from 
these data we can derive upper limits to the total luminos- 
ity of the hot halos. In the special case of the Milky Way 
different pieces of information can be u sed to constrain the 
phys ical properties of the corona (e.g. iFukugita fc Peebles! 
l2006l ). 

Here we consider simple non-rotating spherical models 
of coronae representative of the whim surrounding specific 
disc galaxies, whose gross properties are consistent with the 
observational constraints. The neglect of rotation is a seri- 
ous limitation of this work because rotation must become 
dynamically important sufficiently close to the disc if the 
corona is to feed the disc. Results for non-rotating coronae 
are of interest because the clouds we wish to understand lie 
far from the disc where rotation is expected to be a small 
effect dynamically, and linear analysis of a differentially ro- 
tating coronae is a challenging problem. 

The hot gas is in hydrostatic equilibrium with poly- 
tropic distribution p{r) oc [p(r)]^ , where 7' is the polytropic 
index: the gas temperature profile is 



T(r) 
To 



7' — 1 ^mp 
7' kBTo 



<l>o). 



(14) 



where To = T{ro), <I>o = 'I*('"o) and rg is a reference radius 
- in the following we assume ro — 10 kpc. When 7' > 1 the 
electron number density profile is 



(15) 



where n^o = Mo(''o) and = 0.52n (assuming abundances 
by mass Y = 0.25, X = 0.75). When 7' = 1 the distribution 
is isothermal at temperature To and the electron density 
profile is 



no(r) = no,o exp 



/imp 



$ _ $0) 



(16) 



The gravitational potential "1> is determined by the total 
mass density distribution pdyn, so = 47rG/9dyn. The 

dynamical time is 



tdy 



16Gpdyn 2 ^GMdy„(r) ' 



(17) 



where Mdyn('") = 47r pdyn{r')r'^dr' is the dynamical mass 
of the system and pdyn ~ 3Afdyn('")/47rr^ is the average mass 
density within r. 

We assume here that the total density distributions of 
our model galaxies follow the singular isothermal sphere 
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and that the gravitational potential is 

X. / \ r\ ^B^vir T 

$(r) = 2 In 

/imp \ ro 

where T^ir is the virial temperature, which can be hnked to 
the galaxy's circular speed Vc'. 

± vir — „ , • 



(18) 



(19) 



2kv 



(20) 



When the gas is isothermal (7' = 1) in equilibrium in the 
potential of a singular isothermal sphere, the gas density 
distribution reduces to the simple power law 



no(r) = Hco ( — ) 



(21) 



where a = 2T^ir/To. 

In summary, each galaxy model is characterised by five 
parameters: T^h, 7', To, no,o and the metallicity Z, which 
enters the calculation of the cooling function. 

3.1 Observational constraints 

We use three galaxies as raw models for three categories 
of galaxies with different masses. The first is NGC 5746, a 
very massive spiral with flat rotational speed ~ 310kms~^ 
(corres ponding to a virial tem perature Tvir — 3.4 x 
10^ K: iRand fc Beniamiiil 120081 '). for which we take from 
iRasmussen et al.l (l2009ll the upper limit on the X-ray lu- 
minosity of the corona Lx < 4 x IQ^^ergs"^. The sec- 
ond is the Milky Way, whose corona i s con strained ob- 
servationally as in iFukugita fc PeeblesI l|2006l ). The third 
galaxy is NGC 6503, a low-mass spiral seen at an inclina- 
tion of about 75 degrees, for which Chandra observations 
yield an upper limi t Lx < 1.6 x 10^*ergs~^ to the X- 
ray halo luminosity llStrickland e t al. 2004]). For this galaxy 
we use Vc = 120kms~^ (|Begemani il987l ). which by equa- 
tion (|20[) leads to a virial temperature TVir = 5.1 x 10^ K. In 
both NGC 5746 and in NGC 6503 a component of cxtrapla- 
nar neutral gas has been detected (|Rarid fc Beniamia.2008l : 
ICreisen. Spekkens fc van Moorselll2009l ). 

For each galaxy we consider an isothermal (7' — 1) 
and an adiabatic corona model (7' = 7). We verify that 
the models are consistent with observational constraints by 
computing their X-ray luminosities u sing the Mekal recipe 
l|Liedahl. Osterheld fc Goldsteinll 19951 ) embedded in xspec. 
The 0.3 — 2 keV fluxes are calculated for a range of gas tem- 
peratures (0.1 < T < 1.5 keV) and for different metallici- 
ties (Z). Given fAE{T, Z) the computed speciflc flux, where 
AE = 0.3 — 2keV, the plasma luminosity can be derived as 



Lae{T, Z) = 10" 



dV n^nnfATsiT, Z) ergs~^cm~^ (22) 



where Ue is the electron density and nn is the hydrogen den- 
sity. The factor in front the r.h.s. derives from the xspec nor- 
malization. For isothermal models the temperature in equa- 
tion (I22|) is constant with r and so is fAE{T, Z), so one only 
integrates nl{r), assuming a constant nn/rie ratio. For adi- 
abatic models, the temperature varies with r and /ae(T, Z) 
has to be computed at each r. 

The isothermal model of NGC 5746 has total halo lumi- 
nosity Lx = 3.99 X 10^^ erg s"'^ in the range 5 < r < 40 kpc, 



consistent with the up per limit on the X-ray luminosity 
(jRasmussen et ai1l2009l ) . Our isotherma l model of the Milky 
Way is very similar to "model 2" of iFukugita fc PeeblesI 
(2006). We checked that the obtained luminosity is in agree- 
ment with observational constraints on g alaxies similar to 
the Milky Way, for instance NGC 891 (|Strickland et all 
l2004h . which happens to be the case. For the isother- 
mal model of NGC 6503, we assume nc,o = 1 x 10~^ cm, 
which gives a luminosi ty that matches the upper limit by 
[Strickland et al.1 l|2004h . Similarly, the adiabatic models of 
the three galaxies are constructed so that their X-ray lumi- 
nosities are consistent with observations. 

In contrast to the isothermal case, an adiabatic atmo- 
sphere has a sharp edge and a well defined mass. The lo- 
cation of the corona's edge and thus the total gas mass de- 
pend strongly on the reference gas temperature Tq. We find 
that using the same To (temperature at r = 10 kpc) as the 
isothermal models leads to gas masses which are one or two 
orders of magnitude lower than what is expected for a whim 
corona. Lacking constraints on the temperature profiles, for 
the Milky Way corona we choose To and no,o as the largest 
values (in order to increase the mass of the corona) that do 
not conflict with the luminosity constraints: these values, 
shown in Table[T] give a total whim mass of ~ 1.3 x 10^^ M©, 
in line with the cosmological predictions. For NGC 5746 and 
NGC 6503 we keep the same ratio between To and Tvir as 
for the adiabatic Milky Way model (i.e. same cut-off radius 
~ 350 kpc), maximise ne,o according to the X-ray luminos- 
ity constraints and flnd total masses of the same order (see 
rightmost column of Table [TJ. The density and temperature 
proflles of our models are shown, respectively, in the flrst 
and second rows of panels of Fig. |4] (isothermal models) and 
Fig. [S] (adiabatic models) . The values of the parameters for 
all models are summarised in Table [T] 



4 RESULTS 

4.1 Isothermal coronae 

Consider flrst coronae with isothermal gas distributions 
(7' — 1), which are known to provide reasonable approx- 
imations to the hot atmospheres of elliptical galaxies and 
galaxy clusters. In Fig. [J] we plot the results of applying 
the instability analysis of Section [2.11 to isothermal models 
of NGC 5746 (model N5746i; left-hand column), the Milky 
Way (model MWi; central column) and NGC 6503 (model 
N6503i; right-hand column). The panels in the third row of 
Fig. |4] are the analogues of Fig. |3] and show as functions of 
radius the dynamical, cooling, thermal instability and BV 
timescales. We see that in all cases tev (long-dashed line) is 
shorter than tth (solid line), and the ratio tth/isv is higher in 
galaxies with higher virial temperatures. Hence in all cases 
the growth of thermal perturbations is effectively countered 
by buoyancy. 

A more detailed analysis of the behaviour of the pertur- 
bations can be obtained from the bottom row of panels in 
Fig. |4l which show as a solid line the maximum wavenumber 
ratio (fct /fc)crit a perturbation can have and still grow mono- 
tonically: more tangentially elongated modes [those with 
kt/k > (fct/fc)crit] are overstable. In all cases (fct/fc)crit is 
signiflcantly smaller than unity: this means that it is impos- 
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Table 1. Parameters of the model galaxies. T^i^: virial temperature. 7': polytropic index. Tq: electron temperature at r = lOkpc. n^fi: 
electron number density at r = lOkpc. Z: metallicity. Lx: X-ray luminosity. Afgas: corona gas mass. Isothermal and adiabatic disc galaxy 
models have names ending with "i" and "a", respectively. HM is an isothermal corona model representative of a high-mass elliptical 
galaxy 
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0.51 
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0.82 
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0.03 


0.04 


0.13 


1.01 


N5746a 


3.4 


5/3 


9.7 


0.62 


0.03 


0.27 


12.88 


1.23 


MWa 


1.4 


5/3 


4 


0.68 


0.03 


0.19 


5.58 


1.30 


N6503a 


0.51 


5/3 


1.45 


0.5 


0.03 


0.03 


0.22 


0.96 



sible for such an atmosphere to form blobs (i.e., perturba- 
tions with fcr ~ fct or fct/fc ~ l/\/2) by thermal instability. 
Consistent with the behaviour of the timescales, (fct/fc)crit 
increases for decreasing virial temperature of the system, be- 
ing < 0.007 in model N5746i and 0.06-0.25 in model N6503i. 
Thus, very tangentially elongated disturbances can in prin- 
ciple grow, but we have seen that we do not expect them to 
grow in practice (see Section [2T|. The next-to-bottom row 
of Fig. |4] quantifies ability of thermal conduction to stop 
growth by showing the critical wavelength Acrit as a func- 
tion of radius (solid line). Studies of galaxy clusters indicate 
that a realistic value of the thermal conduc tion suppression 
factor is / ~ 0.01 l|Nipoti fc BinnevI |2004| , and references 
therein), and in this case at a radius r all perturbations 
smaller than ~ O.lr in extent are damped in the Milky- Way 
model MWi. Stabilisation by thermal conduction is propor- 
tionally more effective in higher-temperature systems than 
in lower-temperature systems. For instance, in NGC 5746 at 
all radii perturbations with wave-length smaller than the 
system's radius are damped for / — 0.01, while this is not 
the case in NGC 6503. 

We conclude that if galactic coronae are isothermal, the 
growth of thermal perturbations is effectively suppressed by 
buoyancy and thermal conduction. Suppression is particu- 
larly strong in very massive systems such as NGC 5746, but 
it is very effective also in coronae of Milky- Way-like galax- 
ies and smaller. All these systems are unlikely to form cold 
clouds by thermal instability. 



4.2 Polytropic coronae 

Though isothermality is a well motivated assumption, we 
cannot exclude the possibility that galactic coronae have 
shallower specific entropy profiles, so in this section we in- 
vestigate coronae in which the polytropic relation p{r) oc 
[p(r)]'' holds between pressure and density. As 7' ranges 
from unity to 7 = |, these coronae range from isothermal 
to adiabatic. 

The results for the adiabatic models (N5746a, MWa and 
N6503a in Table 1) are shown in Fig. [5] The cooling and ther- 
mal instability timescales are now not power-law functions 
of radius because density and temperature are no longer 
power-law functions as in the isothermal models (see first 
and second rows of panels) . In an adiabatic atmosphere there 
is no buoyancy to counter the growth of thermal instability, 
so all modes with A > Acrit grow monotonically, independent 



of their wavenumber ratio [tsY — > 00 and (fct/fe)crit — > 00]. 
Thus the thermal stability of the gas is determined by the 
critical wavelength Acrit, which is shown in the bottom row 
of panels in Fig. O If the thermal conduction suppression 
factor is / = 0.01, perturbations are effectively stabilised in 
model N5746a, but not as much in models MWa and N6503a. 
Thus, thermal instabilities cannot develop in the corona of 
NGC 5746, even if it has an adiabatic distribution of gas. 
In the case of the Milky Way, perturbations smaller than 
~ 10 kpc in size are expected to be stabilised by thermal 
conduction for r < 200 kpc. In the case of NGC 6503 pertur- 
bations smaller than 2.5 kpc in size are stable for r < 100 kpc. 
However, if the coronae of lower-mass systems such as the 
Milky Way and NGC 6503 have adiabatic distributions of 
gas, the growth of sufficiently large perturbations by ther- 
mal instability is not excluded at large radii. The rightmost 
panel in the third row of Fig. [5] indicates that at r > 150 kpc 
in model N6503a the perturbation analysis does not apply 
because at these radii fcooi ^ idyn so the corona is not in 
hydrostatic equilibrium; catastrophic cooling is occurring. 

The assumption of perfectly adiabatic atmosphere is 
quite extreme, so it is natural to explore the thermal insta- 
bility of coronae with entropy profiles intermediate between 
the steep ones of isothermal models and the perfectly fiat 
ones of adiabatic models. For this purpose we consider poly- 
tropic models with 1 < 7' < 7. For a polytropic distribution 
with index 7' the BV frequency is 



2 _ g fJ-rrip 



T ks 



where 



/i(7,7') 



1/2 



(23) 



(24) 



Thus we obtain the following expression for the maximum 
ratio of wavenumbers that is consistent with growing mode: 



(-) 

\ k J crit 



7- 



1 



nA{T)A{T) 



27/1(7,7') (^mpfeB)i/2Ti/2g- (2^) 

We focus on Milky- Way-like polytropic models of index 7', 
with Tvir = 1.4 X 10^ K, Z = O.OS^'q, and values of nc,o and 
To such that they have the same gas mass within 500 kpc 
and the same value of the specific entropy at the radius 
Thaif containing half of the gas mass as the adiabatic model 
MWa. In Fig. |S]we plot as a solid line (fct/fc)crit as a func- 
tion of 7' at r = Thaif, which is in the range 140 — 170 kpc 
for these models. Only perturbations with kt/k < (fct/fc)crit 
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Figure 4. Prom top to bottom: electron density, temperature, timescales, perturbation wavelength, and tangential-to-total perturbation 
wavenumber ratio as functions of radius for the isothermal models N5746i (left-hand column), MWi (central column) and N6503i (right- 
hand column). The timescales panels plot the dynamical (short dashed), BV (long dashed), cooling (dotted) and thermal instability (solid) 
timescales. In the wavelength panels the solid curves represent the critical wavelength A^rit for stabilisation by thermal conduction: only 
perturbations with wavelengths in the shaded areas grow (/ is the factor by which thermal conduction is suppressed; see equation |3)l . In 
the wavenumber ratio panels the solid curves represent the critical wavenumber ratio (fct/fc)crit for instability: only perturbations with 
ratio in the shaded areas are unstable. 



(shaded area in the diagram) are unstable. We see that as 
soon as 7' deviates from 7, (fct/fc)crit drops from infinity to 
values smaller than unity, and (fet/fc)crit ^ 0.2 for 7' < 1.55. 
This means that thermal instability is efficiently countered 
by buoyancy even in polytropic distributions with quite shal- 
low entropy profiles. We conclude that unless the Milky Way 



corona is very close to adiabatic, thermal perturbations can- 
not grow. 

4.3 Time evolution of polytropic coronae 

In the absence of strong observational constraints, we try 
here to address on theoretical grounds the question of 
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Figure 5. Same as Fig. |4] but for the adiabatic models N5746a (left-hand column), MWa (central column) and N6503a (right-hand 
column). In these cases Ibv and (fct/A;)crit are infinite at all radii, so we do not plot diagrams of kt/k. 



whether the coronae of disk galaxies are close adiabatic 
by exploring how a corona with a shallow entropy profile 
evolves in time in the presence of cooling. If a steep outward- 
increasing entropy profile were produced quickly, we could 
discard the hypothesis that coronae have very shallow en- 
tropy profiles. For this purpose, we consider the extreme case 
of a corona that starts from a perfectly flat entropy proflle 
('y' = 'y polytrope) and the more general case of coronae 
with shallow outward increasing entropy profiles (7' < 7 
polytropes) . 

The left panels in Fig.[7]show how an initially adiabatic 
model of the corona of the Milky Way evolves over 1 Gyr. 
Although the central density rises steadily (top left), the 
temperature profile changes significantly only at the outer 
edge of the corona (bottom left). Here the temperature was 
initially 3 x 10^ K, and, as one moves outwards, the in- 
crease in the efficiency of cooling with falling temperature 
outweighs the reduction in efficiency with decreasing den- 
sity. Consequently, the initially flat profile of speciflc en- 



tropy a starts to fall with radius. In this region of outwards- 
decreasing speciflc entropy, the corona is convectively un- 
stable and clouds of cool gas will fall inwards. Further in, 
the speciflc entropy proflle remains extremely flat as entropy 
is radiated away, and regions of depressed temperature will 
not be stabilised by buoyancy. Thus a non-linear computa- 
tion conflrms that in the extreme case of an adiabatic at- 
mosphere perturbations large enough not to be damped by 
thermal conduction can form cool clouds. 

It is unlikely that all the corona's gas will have exactly 
the same speciflc entropy - initially some gas is bound to 
have more specific entropy than other gas, and the galaxy's 
gravitational field will sort the gas by speciflc entropy such 
that low-entropy gas lies inside higher-entropy gas. We now 
address the question of how wide a spread in speciflc en- 
tropy is required to prevent cooling producing a region of 
outwards-decreasing speciflc entropy in which clouds could 
form. A simple way of answering this question is to consider 
the cooling of coronae that are initially polytropes of index 
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Figure 7. Evolution over 1 Gyr of model coronae for the Milky Way. Left column: an initially adiabatic model (7' = 1.6667). Right 
column: a model with 7' = 1.55. The full curves in the top panels show the increase in the central electron density over 1 Gyr, while the 
dotted curve in the top right panel shows the initial density profile of the adiabatic model of the left column. In the central and bottom 
panels the left scale gives a = pp^"' (in units of kcVcm^ ) and the right scale gives logT. In the central panels the independent variable 
is radius and broken curves show temperature (the dotted curve shows the initial condition). In the bottom panels the independent 
variable is mass shell and initial and final conditions are shown by full and dashed curves respectively. The coronae were evolved by the 
technique described in Kaiser & Binney (2003) with metallicity Z = O.OSZq. 



7' < 7, and we focus here on the Milky- Way-hke models of 
Fig.H 

The right panels of Fig. [7] show the evolution of the 
Milky- Way corona that starts from the polytrope with 7' = 
1.55. The full curve labelled "ct" in the bottom right panel 
shows this corona initially has a significant entropy gradi- 
ent in the outer region, where an entropy depression formed 
in the adiabatic case. The dotted curve for a in the bot- 
tom right panel shows that after 1 Gyr, the entropy gradient 
has been reversed, and there is again a region of outwards- 
decreasing specific entropy. For polytropic indices 7' < 1.4 a 



cooling catastrophe develops at the centre before there has 
been significant cooling at the periphery, and for 7' = 1.4 a 
point of inflection develops in cr(M) at M = 12.2 x 10^° Mq 
around the time (~ 0.5 Gyr) when cooling becomes catas- 
trophically fast at the centre. 

The central panels of Fig. [7] show that thermal instabil- 
ity occurs only at radii in excess of 100 kpc. This is the case 
because the steep portion in the cooling curve, which drives 
the instability, occurs at temperatures that are significantly 
below the central temperature of the corona - such temper- 
atures are reached only close to the edge of the corona be- 
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Figure 6. Critical wavenumber ratio (fet/fc)crit as a function of 
the polytropic index 7' at the half-mass radius r^aif for Milky- 
Way-hke polytropic models with Tvir = 1.4 X 10^ K and Z = 
O.OSZq. 



cause when the circular speed is assumed constant (as here) 
they imply a steep density profile. In fact, the formation of a 
region of outwards-decreasing specific entropy depends sen- 
sitively on the shape of the cooling function. When A(r) is 
that given by Sutherland & Dopita (1993) for zero metallic- 
ity, no such region forms because A(r) does not fall steeply 
enough with increasing temperature. At zero metallicity an 
initially adiabatic corona develops an outwards-increasing 
entropy gradient on account of the increase in the radiation 
rate with increasing density. 

Once entropy is either independent of radius or 
outwards-decreasing, there is nothing to inhibit bulk radial 
motion of gas. In an adiabatic atmosphere such motions need 
to be externally excited, for example by tidal fields and the 
motion of dark-matter substructures and satellite galaxies 
as they orbit through the host gravitational potential. In a 
convectively unstable atmosphere, such motions will arise 
spontaneously. Radial bulk motions transfer heat between 
different radial shells of a corona, and, as in the theory of 
stellar structure, it is natural to assume that convective en- 
ergy transport will maintain an adiabatic entropy profile in 
any region in which entropy would otherwise be outwards- 
decreasing. Hence we modified the code that produced Fig. [7] 
to include convection as follows. After each timestep, the 
code looks for regions in which entropy is outwards decreas- 
ing and sets the specific entropy Si = ln(a)+constant of each 
mass shell in this region to a common value s. The search 
for regions of outwards-decreasing entropy, and subsequent 
entropy adjustment is then repeated until no such region is 
found. The entropy s is reached as a result of the ith mass 
shell receiving an amount of heat dqi from the other shells, 
thus changing its specific entropy by 




Figure 8. The evolution of the initial polytrope with 7' = 1.55 
when convection is assumed to prevent entropy decreasing out- 
wards. Top panel as a function of mass; lower panel as a function 
of radius. 
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where Ti and Mi are the temperature and mass of the ith 
shell. By conservation of energy, ^ , dq^ = 0, so after mul- 
tiplying (1261) through by TiMi and summing over i we have 

Fig. |8] shows the evolution of the initially polytropic 
corona plotted in the right panels of Fig. [7] when this algo- 
rithm is used to simulate the effects of convection. One sees 
that after a Gyr specific entropy is constant through nearly 
half the mass. Thus metal-enhanced cooling combined with 
convection has considerable potential for driving the outer 
half of a polytropic corona to an adiabatic configuration, 
provided that the initial entropy profile is shallow enough. 
However, the fact that polytropic models are characterised 
by a finite boundary at which the temperature drops to zero 
is crucial for these systems to evolve towards a configuration 
with fiat entropy profile in the outer regions. We do not ex- 
pect that these conditions are met in real galaxies, which 
are surrounded by the intergalactic medium. 



5 CONCLUSIONS 

In this paper we examined the proposal that HVCs of the 
Milky Way and other disc galaxies form by condensation of 
the WHIM through thermal instability. The linearised equa- 
tions for the evolution of thermal perturbations within a 
cooling, conducting, non-rotating spherical corona show that 
unless specific entropy is a shallow function of radius (a 
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nearly adiabatic corona), regions in which the temperature 
is initially anomalously low do not experience the classical 
Field instability. In fact they are simply points at which an 
internal gravity wave is momentarily in the ascendant, and 
after half a Brunt- Vaisala period the material in this region 
will sink to below its equilibrium radius and will then be an 
anomalously warm region. Consequently, there is no thermal 
instability from which clouds can form in a corona in which 
there is a significant gradient in the specific entropy. 

Clouds can form in a strictly adiabatic corona because 
in this case buoyancy vanishes and the increase in the spe- 
cific cooling function A(r) with decreasing temperature 
around T ~ 3 x 10^ K can cause the cooling rate to in- 
crease outwards notwithstanding the steady decline in gas 
density with radius. This increase in the cooling rate can 
cause the specific entropy to become outwards-decreasing 
and thus the corona to become convectively unstable. In 
these circumstances clouds of rapidly cooling material will 
tumble inwards if thermal conductivity is too weak to carry 
enough heat from the warm interior to the rapidly cooling 
periphery to offset radiative cooling. 

An examination of adiabatic coronae for three galax- 
ies, NGC 5746, the Milky Way, and NGC 6503 shows that 
thermal conductivity that is suppressed to 1% of Spitzer's 
fiducial value would eliminate thermal instability entirely 
in NGC 5746 and largely in the Milky Way, where pertur- 
bations smaller than lOkpc in extent could grow only at 
r > 200 kpc. The efficiency with which thermal conductivity 
can suppress thermal instability is greatest in the hottest 
coronae because the Spitzer's conductivity is a strongly in- 
creasing function of temperature. So in a low-mass galaxy 
such as NGC 6503, perturbations 10 kpc in extent could grow 
through most of the corona, although a perturbation only a 
kpc in extent would be stabilised almost everywhere. 

It is implausible that all the gas in a corona initially has 
exactly the same specific entropy. The gas will inevitably 
be spread over a range of specific entropies, and within 
a few dynamical times convection will sort the corona by 
specific entropy so that entropy is outwards-increasing. A 
convenient way of simulating the consequences of such a 
spread in entropy is to consider coronae that are initially 
polytropes with index 7' < 7. In such coronae there is a sig- 
nificant radial gradient in specific entropy near the corona's 
edge, where the temperature is low and thermal instability 
is most likely. Non-linear simulations of the cooling of ini- 
tially polytropic coronae for the Milky Way show that when 
the polytropic index 7' is smaller than 1.4, cooling steepens 
the initial outwards-increasing entropy gradient rather than 
reversing it and this ensures that anomalously cool regions 
are buoyantly stabilised against thermal instability. 

Our conclusion that a flat entropy proflle is a necessary 
condition for the growth of thermal instability in galactic 
coronae is consistent with the re cent results of SPH simula- 
tions bv lKaufmann et al.1 l|2009l '). who find that cold clouds 
form from hot gas only in systems with flat entropy pro- 
flies. However, we find that except in very low-mass sys- 
tems, thermal conductivity (which was not considered by 
[Kaufmann ct al. 2009), is capable of suppressing thermal 
instability at most radii, even when the conductivity is 
strongly suppressed. 

In light of these results, is it likely that the clouds of H I 
that are actually observed have formed by cooling of coro- 



nal gas? A major issue is that observed H I clouds, such 
as Comple x C and the filament of NGC 891 lie at radii 
~ 10 kpc (|Qosterloo. Fraternah fc Sancisi ] l2007l ). Clouds 
are not found as far out as the bounding radii of a nearly adi- 
abatic atmospheres that would contain the missing baryons, 
while theory indicates that it is precisely near these radii 
that clouds may form as a result of the coronal temperature 
falling through T ~ 3 x 10* K, where the cooling function is 
sufficiently steep. One might speculate that initially clouds 
have densities and masses that lie below the thresholds for 
detection, and it is only after they have fallen in to much 
smaller radii, where the confining pressure is higher, and ac- 
creted coronal gas or merged with other clouds, that they 
become detectable. In this case deeper surveys should find 
clouds out to in excess of 100 kpc. However, a difficulty with 
this proposal is that small clouds are thermally unstable only 
if the corona's thermal conductivity is heavily suppressed - 
in fact to less than a percent of Spitzer's value. 

Thermal instability only occurs if the coronal temper- 
ature falls through T ~ 3 x 10* K. In the polytropic coro- 
nae discussed here this condition is satisfied because T van- 
ishes at a finite radius. However, this feature of a polytropic 
corona is artificial: in reality the corona will at some point 
merge into the general intergalactic medium, which is likely 
to have high entropy and a temperature which could well be 
above 5 x 10* K. 

Except in Section r4.3l our analysis is linear in nature and 
does not apply close to the star-forming disc since this region 
is violently disturbed by expanding supershells. However, 
linear analysis should be a reasonable approximation several 
kiloparsecs from the disc where HVCs are observed. 

A significant limitation of our work is th at it excludes 
both magn etic field s and r otation of the corona. [Loewensteinl 
(1990) and iBalbud (|l99lD showed that magnetic fields can 
in principle have a big impact on the stability of cooling 
fiows. However, the relevance of the linear calculations pre- 
sented by these authors to real cooling flows and coronae is 
doubtful because in order to perform the analysis one has 
to assume an unperturbed magnetic-field configuration that 
is consistent with the unperturbed cooling flow - this field 
turns out to be radial. The modes obtained obviously reflect 
the particular stress tensor of this fleld, and are valid only 
for perturbations that do not signiflcantly change the fleld 
geometry because (5B ^ B. In cooling flows magnetic pres- 
sure is much smaller than thermal pressure and probably on 
the order of the turbulent pr essure put^^j,. Consequently the 
results of lLoewensteinI l\199dl ) are valid only when the veloci- 
ties associated with modes are small compared to «turb- But 
we are precisely interested in the turbulent motions them- 
selves, so this condition will be seriously violated. Moreover, 
given that the magnetic pressure is not larger than the tur- 
bulent pressure, the magnetic fleld will be tangled by the 
turbulence rather than have the special geometry required 
by the assumption that a steady state is maintained in the 
presence of a radial accretion flow. It seems likely that in 
these circumstances the safest procedure is simply to assume 
that the effect of the (tangled) magnetic fleld is conflned to 
augmenting the isotropic pressure by ~ /2fio- In this case 
our results are applicable. 

Rotation, by contrast, may have a signiflcant impact on 
thermal stability. Indeed, if clouds that condense from the 
corona are to feed the growth of a centrifugally supported 
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gas disc, rotation must be dynamically important near the 
disc. However, clouds are observed far from the disc, where 
rotation will be less important, and the question is whether 
a level of rotation that is dynamically small (in the sense 
that the centrifugal acceleration is much smaller than the 
gravitational acceleration) can introduce thermal instabil- 
ity by obstructing internal gravity waves. We reserve this 
difficult question for a later paper. 

We conclude that, if galactic haloes are well represented 
by non-rotating, quasi-hydrostatic atmospheres, thermal in- 
stability is not a viable mechanism for the formation of cold 
clouds in disc galaxies as massive as, or more massive than, 
the Milky Way. It may just work in low-mass systems if 
their (essentially unconstrained) entropy profiles are almost 
perfectly fiat. Our best guess is that the observed popula- 
tion of HVCs around the Milky Way and other large nearby 
galaxies represent gas that has been stripped from satel- 
lites. This conclusion leaves open the question of how gas 
moves from the whim t o the disc, as it must if sta r forma- 
tion is to be sustained. iFraternali fc BinnevI (|2008l ) present 
circumstantial evidence that the key to this transfer is the 
galactic fountain, which constantly cycles disc gas through 
the bottom of the corona. To understand this transfer we 
need models of the corona that include significant centrifu- 
gal support near the disc. 
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